*********************************************************************************
*	Polarization Between the Rich and the Poor?									*
*	The Dynamics and Structure of Redistributive Preferences in a Comparative	*
*	Perspective	(supplementary materials)										*
*																				*
*	Xavier Romero-Vidal, University of Cambridge								*
* 	Steven M. Van Hauwaert, University of Surrey & Radboud University Nijmegen	*
*********************************************************************************

encode var, gen(var2)
xtset var2 year

* D.1 and D.2
alpha std_fr_all std_de_all std_gb_all std_nl_all std_no_all std_se_all std_ch_all, detail item

* E.1 and E.2 
factor std_fr_all std_de_all std_gb_all std_nl_all std_no_all std_se_all std_ch_all

*F: 
alpha std_fr_poor std_fr_mid std_fr_rich, item detail
alpha std_de_poor std_de_mid std_de_rich, item detail
alpha std_gb_poor std_gb_mid std_gb_rich, item detail
alpha std_nl_poor std_nl_mid std_nl_rich, item detail
alpha std_no_poor std_no_mid std_no_rich, item detail
alpha std_se_poor std_se_mid std_se_rich, item detail
alpha std_ch_poor std_ch_mid std_ch_rich, item detail

*G (in a dataset where countries are columns): 
factor std_fr_poor std_de_poor std_gb_poor std_nl_poor std_no_poor std_se_poor std_ch_poor
factor std_fr_mid std_de_mid std_gb_mid std_nl_mid std_no_mid std_se_mid std_ch_mid
factor std_fr_rich std_de_rich std_gb_rich std_nl_rich std_no_rich std_se_rich std_ch_rich

*H
anova mood year if class != 0 // .145
anova mood year if country==1 & class != 0 // France (.556)
anova mood year if country==2 & class != 0 // Germany (.372)
anova mood year if country==3 & class != 0 // Great Britain (.634)
anova mood year if country==4 & class != 0 // Netherlands (.593)
anova mood year if country==5 & class != 0 // Norway (.271)
anova mood year if country==6 & class != 0 // Sweden (.439)
anova mood year if country==7 & class != 0 // Switzerland (.241)

* Variance associated to income group

anova mood class if class != 0 // .277
anova mood class if country==1 & class != 0 // France (.285)
anova mood class if country==2 & class != 0 // Germany (.312)
anova mood class if country==3 & class != 0 // Great Britain (.288)
anova mood class if country==4 & class != 0 // Netherlands (.231)
anova mood class if country==5 & class != 0 // Norway (.551)
anova mood class if country==6 & class != 0 // Sweden (.503)
anova mood class if country==7 & class != 0 // Switzerland (.271)

*I.1
* FE model
xtreg std_mood_mid l.std_mood_mid std_mood_poor std_mood_rich, fe
* RE model
xtreg std_mood_mid l.std_mood_mid std_mood_poor std_mood_rich, re
* GLS model
xtgls std_mood_mid l.std_mood_mid std_mood_poor std_mood_rich, corr(ar1)
* PRAIS-Winsten model
prais std_mood_mid l.std_mood_mid std_mood_poor std_mood_rich, corc
* PCSE model
xtpcse std_mood_mid l.std_mood_mid std_mood_poor std_mood_rich, corr(psar1)

*Model 7:
* FE model
xtreg d.std_mood_mid l.std_mood_mid d.std_mood_poor d.std_mood_rich, fe
* RE model
xtreg d.std_mood_mid l.std_mood_mid d.std_mood_poor d.std_mood_rich, re
* GLS model
xtgls d.std_mood_mid l.std_mood_mid d.std_mood_poor d.std_mood_rich, corr(ar1)
* PRAIS-Winsten model
prais d.std_mood_mid l.std_mood_mid d.std_mood_poor d.std_mood_rich, corc
* PCSE model
xtpcse d.std_mood_mid l.std_mood_mid d.std_mood_poor d.std_mood_rich, corr(psar1)

***************************************


*Switzerland
xtset var2 year

twoway (tsline value if var2==1, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==2, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==3, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==4, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==5, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==6, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==8, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==9, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==10, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==11, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==12, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==13, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==14, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==15, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==16, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==17, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==18, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==19, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==20, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==21, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==22, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==23, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==24, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==25, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==26, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==27, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==28, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==29, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==30, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==31, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==32, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==33, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==34, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==23, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==36, lcolor(gs12) lwidth(vvvthin))  ///
(tsline value if var2==7, lcolor(black) lwidth(medium)), ///
scheme(plottig) ytitle(Redistributive preferences, size(small)) xlabel(1993(4)2017, angle(45) labsize(vsmall)) graphregion(color(white)) bgcolor(gs16) legend(off)

*Germany
xtset var2 year

twoway (tsline value if var2==1, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==2, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==3, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==4, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==5, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==6, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==8, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==9, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==10, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==11, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==12, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==13, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==14, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==15, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==16, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==17, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==18, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==19, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==20, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==21, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==22, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==23, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==24, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==25, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==26, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==27, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==28, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==29, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==30, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==31, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==32, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==33, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==34, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==35, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==36, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==37, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==38, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==39, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==40, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==41, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==42, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==43, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==44, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==45, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==46, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==47, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==48, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==49, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==50, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==51, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==52, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==53, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==24, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==55, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==56, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==57, lcolor(gs12) lwidth(vvvthin)) ///
														 (tsline value if var2==59, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==60, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==61, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==62, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==63, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==64, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==65, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==66, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==67, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==68, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==69, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==70, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==71, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==72, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==73, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==74, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==75, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==58, lcolor(black) lwidth(medium)), ///
scheme(plottig) ytitle(Redistributive preferences, size(small)) xlabel(1980(5)2015, angle(45) labsize(vsmall)) graphregion(color(white)) bgcolor(gs16) legend(off)


*France
encode var, gen(var2)
xtset var2 year
twoway (tsline value if var2==1, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==2, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==3, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==4, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==5, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==6, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==8, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==9, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==10, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==11, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==12, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==13, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==14, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==15, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==16, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==17, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==18, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==19, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==20, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==21, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==22, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==23, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==24, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==25, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==26, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==27, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==28, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==29, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==30, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==31, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==32, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==33, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==34, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==35, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==36, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==37, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==38, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==39, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==40, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==41, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==42, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==43, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==44, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==45, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==46, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==47, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==48, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==49, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==50, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==51, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==52, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==53, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==24, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==55, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==56, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==57, lcolor(gs12) lwidth(vvvthin)) ///
														 (tsline value if var2==59, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==60, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==61, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==62, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==63, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==64, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==65, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==66, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==67, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==68, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==69, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==70, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==71, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==72, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==73, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==74, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==75, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==76, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==77, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==78, lcolor(gs12) lwidth(vvvthin))  ///
(tsline value if var2==58, lcolor(black) lwidth(medium)), ///
scheme(plottig) ytitle(Redistributive preferences, size(small)) xlabel(1980(5)2015, angle(45) labsize(vsmall)) graphregion(color(white)) bgcolor(gs16) legend(off)

*GB
twoway (tsline value if var2==1, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==2, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==3, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==4, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==5, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==6, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==8, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==9, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==10, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==11, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==12, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==13, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==14, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==15, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==16, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==17, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==18, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==19, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==20, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==21, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==22, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==23, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==24, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==25, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==26, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==27, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==28, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==29, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==30, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==31, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==32, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==33, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==34, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==35, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==36, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==37, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==38, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==39, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==40, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==41, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==42, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==43, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==44, lcolor(gs12) lwidth(vvvthin)) 														  ///
(tsline value if var2==46, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==47, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==48, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==49, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==50, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==51, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==52, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==53, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==24, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==55, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==56, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==57, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==58, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==59, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==60, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==61, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==62, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==63, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==45, lcolor(black) lwidth(medium)), ///
scheme(plottig) ytitle(Redistributive preferences, size(small)) xlabel(1980(5)2015, angle(45) labsize(vsmall)) graphregion(color(white)) bgcolor(gs16) legend(off)

*Netherlands
twoway (tsline value if var2==1, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==2, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==3, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==4, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==5, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==6, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==8, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==9, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==10, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==11, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==12, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==13, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==14, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==15, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==16, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==17, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==18, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==19, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==20, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==21, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==22, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==23, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==24, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==25, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==26, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==27, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==28, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==29, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==30, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==31, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==32, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==33, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==34, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==35, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==36, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==37, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==38, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==39, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==40, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==41, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==42, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==43, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==44, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==45, lcolor(gs12) lwidth(vvvthin)) ///														  ///
(tsline value if var2==46, lcolor(gs12) lwidth(vvvthin)) 														  ///
(tsline value if var2==48, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==49, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==50, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==51, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==52, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==53, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==24, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==55, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==56, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==57, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==58, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==59, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==60, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==61, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==62, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==63, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==47, lcolor(black) lwidth(medium)), ///
scheme(plottig) ytitle(Redistributive preferences, size(small)) xlabel(1980(5)2015, angle(45) labsize(vsmall)) graphregion(color(white)) bgcolor(gs16) legend(off)

*Sweden
twoway (tsline value if var2==1, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==2, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==3, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==4, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==5, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==6, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==7, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==10, lcolor(gs12) lwidth(vvvthin)) /// 
(tsline value if var2==11, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==12, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==13, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==14, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==15, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==16, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==17, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==18, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==19, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==20, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==21, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==22, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==23, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==24, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==25, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==26, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==27, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==28, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==29, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==30, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==31, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==32, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==33, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==34, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==35, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==36, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==37, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==38, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==39, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==40, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==41, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==42, lcolor(gs12) lwidth(vvvthin)) ///
														(tsline value if var2==44, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==45, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==46, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==47, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==48, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==49, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==50, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==51, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==52, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==53, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==54, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==55, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==56, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==57, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==58, lcolor(gs12) lwidth(vvvthin)) ///		
(tsline value if var2==59, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==60, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==61, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==62, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==63, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==64, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==65, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==66, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==67, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==68, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==69, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==70, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==71, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==72, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==73, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==74, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==43, lcolor(black) lwidth(medium)), ///
scheme(plottig) ytitle(Redistributive preferences, size(small)) xlabel(1980(5)2015, angle(45) labsize(vsmall)) graphregion(color(white)) bgcolor(gs16) legend(off)

*Denmark
twoway (tsline value if var2==1, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==2, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==3, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==4, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==5, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==6, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==7, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==10, lcolor(gs12) lwidth(vvvthin)) /// 
(tsline value if var2==11, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==12, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==13, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==14, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==15, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==16, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==17, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==18, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==19, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==20, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==21, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==22, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==23, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==24, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==25, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==26, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==27, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==28, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==29, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==30, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==31, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==32, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==33, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==34, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==35, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==36, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==37, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==38, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==39, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==40, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==41, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==42, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==33, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==44, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==45, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==46, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==47, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==48, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==49, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==50, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==51, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==52, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==53, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==54, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==55, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==56, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==57, lcolor(gs12) lwidth(vvvthin))														 ///
(tsline value if var2==59, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==60, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==61, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==62, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==63, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==64, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==65, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==66, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==67, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==68, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==69, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==70, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==71, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==72, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==73, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==74, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==58, lcolor(black) lwidth(medium)), ///
scheme(plottig) ytitle(Redistributive preferences, size(small)) xlabel(1980(5)2015, angle(45) labsize(vsmall)) graphregion(color(white)) bgcolor(gs16) legend(off)

*Norway
xtset var2 year
twoway (tsline value if var2==1, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==2, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==3, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==4, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==5, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==6, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==7, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==8, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==9, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==10, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==11, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==13, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==14, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==15, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==15, lcolor(gs12) lwidth(vvvthin)) (tsline value if var2==16, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==17, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==18, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==19, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==20, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==21, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==22, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==23, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==24, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==25, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==26, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==27, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==28, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==29, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==30, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==31, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==32, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==33, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==34, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==35, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==36, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==37, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==38, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==39, lcolor(gs12) lwidth(vvvthin))(tsline value if var2==40, lcolor(gs12) lwidth(vvvthin)) ///
(tsline value if var2==12, lcolor(black) lwidth(medium)), ///
scheme(plottig) ytitle(Redistributive preferences, size(small)) xlabel(1981(4)2016, angle(45) labsize(vsmall)) graphregion(color(white)) bgcolor(gs16) legend(off)


*Panelled figure
graph combine FR.gph DE.gph GB.gph NL.gph NO.gph SE.gph CH.gph, ///
scheme(plottig) 
